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Abstract 



Rate-independent systems allow for solutions with jumps that need additional 
modeling. Here we suggest a formulation that arises as limit of viscous regulariza- 
tion of the solutions in the extended state space. Hence, our parametrized metric 
solutions of a rate-independent system are absolutely continuous mappings from 
^ ' a parameter interval into the extended state space. Jumps appear as generalized 

^^ ■ gradient flows during which the time is constant. The closely related notion of BV 

[■^^ . solutions is developed afterwards. Our approach is based on the abstract theory of 

^^ I generalized gradient flows in metric spaces, and comparison with other notions of 

C^^ ' solutions is given. 

o- 
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X ■ 1 Introduction 



This paper is concerned with the analysis of different solution notions for rate- independent 
evolutionary systems. The latter arise in a very broad class of mechanical problems, 
usually in connection with hysteretic behavior. With no claim at completeness, we may 
mention for instance elastoplasticity, damage, the quasistatic evolution of fractures, shape 
memory alloys, delamination and ferromagnetism, referring to |Mie05] for a survey of the 
modeling of rate-independent phenomena. 

Because of their relevance in applications, the analysis of these systems has attracted 
some attention over the last decade, also in connection with the issue of their proper 
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formulation. In fact, in several situations rate- independent problems may be recast in the 
form of a doubly nonlinear evolution equation involving two energy junctionals, namely 

dq'J{.{q{t),q{t)) + dqE{t,q{t))3Q in Q' for a.a. t G (0, T) , (1.1) 

where Q is a separable Banach space, 3^ : Q x Q — > [0, oo] a dissipation functional and 
£ : [0, T] X Q — >• (— oo, oo] an energy potential, dq and dq denoting their sub differential with 
respect to the second variable. Rate-independence is rendered through 1-homogeneity of 
the functional 51 with respect to its second variable. Indeed, assuming that 3l(g, 7f ) = 
73i(g, f ) for all 7 > and (g, f) E Q x Q, one has that equation f 1 1.11) is invariant for 
time-rescalings. This captures the main feature of this kind of processes, which are driven 
by an external loading set on a time scale much slower than the time scale intrinsic to the 
system, but still fast enough to prevent equilibrium. Typically, this quasistatic behavior 
originates in the limit of systems with a viscous, rate-dependent dissipation. 

The formulation of rate-independent problems in terms of the sub differential inclusion 
(11. ip has been thoroughly analyzed in [MiT04j . in the case of a reflexive Banach space. 
Existence of solutions to the Cauchy problem for (11.11) is proved through approximation 
by time discretization and solution of incremental minimization problems. However, in 
many applications the energy £(t, ■) is neither smooth nor convex, and the state space 
Q is often neither reflexive nor the dual of a separable Banach space (see e.g. [KrZ07j 
for energies having linear growth at infinity). Furthermore, Q may even lack a linear 
structure (in these cases we will denote it by the calligraphic letter Q), like for finite- 
strain elastoplasticity |Mie03t IMaMOSj or for quasistatic evolution of fractures |DFT05] . 

In such situations, the differential formulation (11.11) cannot be used. In |MiT99l 
IMTL02J . the concept of energetic solution for general rate-independent energetic systems 
(Q, £, T)) has been introduced, by replacing the infinitesimal metric CR of the subdifferential 
formulation (II. ip by a global dissipation distance 2) : Q x Q — > [0, oo]. This formulation, 
see Section 15.11 is derivative-free, and thus applies to solutions with jumps and can be 
used in very general frameworks, like, for example, in a topological space Q, with £ and 2) 
lower semi-continuous only, see |MaM05( IMieOSl [FrMOGj . Energetic solutions are very flex- 
ible and allow for a quite general existence theory; however, the global stability condition, 
asking that q(t) globally minimizes the map q\-^ £(t, q)-\-T){q(t),q), implies that solutions 
jump earlier as physically expected, since they are forced to leave a locally stable state, see 
e.g. |Mie03t Ex. 6.1] or [KMZ07t Ex. 6.3], and Example 17.11 below. Moreover, existence of 
energetic solutions is proved via time discretization and incremental global minimization, 
but, as discussed in |Mie03l Sec. 6], local minimization would be more appropriate both 
from the perspective of modeling and of numerical algorithms. 

In response to these issues, in |EfM06] a vanishing viscosity approach was proposed 
to derive new solution types for rate-independent systems (Q, "R, £). There, Q is assumed 
to be a finite- dimensional Hilbert space Q and £ e C^([0,T] x Q). The natural viscous 
approximation of (11.11) is obtained by adding a quadratic term to the dissipation potential, 
viz. 3ie(g, t>) = 'Ji{q,v) + §11^11^, and leads to the doubly nonlinear equation 

Eq{t)+dgJl{q{t),q{t)) + dqE{t,q{t))3 for a.a. t G (0, T) . (1.2) 



Using dim Q < oo, the existence of solutions g^ G H^([0,T];(5) is obvious and passing 
to the hmit e \ in fll.2p leads to new solutions and to a finer description of the 
jumps, which occur later than for energetic solutions. The key idea (see Section [2]) is 
that the limiting solution at jumps shall follow a path which somehow keeps track of 
the viscous approximation. To exploit this additional information, one has to go over 
to an extended state space: reparametrizing the approximating viscous solutions q^ of 
(11. 2p by their arclength r^, and introducing the rescalings t^ = t^^ and q^ = qe °te, one 
studies the limiting behavior of the sequence {(te, qe)}e as e | 0. Hence, in [EfMOGj it was 
proved that (up to a subsequence), {(te,^^)}^ converges to a pair (t,q), whose evolution 
encompasses both dry friction effects and, when the system jumps, the infiuence of rate- 
dependent dissipation. In fact, the jump path may be completely described by a gradient 
fiow equation, which leads to this interpretation: jumps are fast (with respect to the slow 
external time scale) transitions between two metastable states, during which the system 
switches to a viscous regime. Furthermore, solutions of the limiting rate-independent 
problem can be constructed by means of a time-discretization scheme featuring local, 
rather than global, minimization. 

This paper provides the first step of the generalization of these ideas to the much more 
general metric framework using the concept of curves of maximal slope, which dates back 
to the pioneering paper |DGMT80] . We also refer to the recent monograph [AGS05J . the 
references therein, and to [RMSOSj . The general setup starts with a 

complete metric space (X, d) 

and introduces the metric velocity 

h\0 h h\0 h 

which is defined a.e. along an absolutely continuous curve q : [0, T] -^ X. 

This replaces the norm of the derivative q' in the smooth setting, and, in the same 
way, the norm of the (Gateaux)-derivative or the sub differential of a functional 3" : X ^ 
(— oo, oo] is replaced by the local slope of 3" in g G dom(\E'), which is defined by 

\d,n<i) :=iimsup^^^4r^r^' (1-4) 

where (■)^ denotes the positive part. With these concepts, the viscous problem (11.21) has 
the equivalent metric formulation 



d 
dt 



E,{t,q{t))-dMt,q{t)) < -iWm + IWm) - ^ {{\d,E,\ (t,g(t)) - l)+)', (1.5) 



for a.a. t G (0,T), see Section ISTTl for further details. It was proved in [RMSOSl Thm. 3.5] 
that, under suitable assumptions on £, for every e > the related Cauchy problem has 
at least one solution q^ G AC([0,T]; X). 

Following the approach of |EfM06] . for every e > we now consider the (arclength) 
rescalings (te,^^) associated with q^, which in turn fulfill a rescaled version of (11.51) . cf. 



f l3.12p . Under suitable assumptions, in Theorem 13.81 we shall show that, up to a sub- 
sequence, {Cts,(le)} converges to a limit curve (t,g) G AC([0, S']; [0,T] x X) such that 



i : [0, S] -^ [0,T] is nondecreasing, 
?(s) + \q\{s) > for a.a. s E [0, S], 

?(s)>0 =^ |9/|(f(s),g(s))<l, 
|g'|(s)>0 ^ \d,^\{t{s),q{s))>l, 



(1.6a) 
for a.a. s G [0, S] , (1.6b) 



and the energy identity 



Ae(f(,),g(,))-9,£(f(s),g(s))?(s) = (D,£(f(.),g(.)),g'(.)) ^^ ^^^ 

= -\q\{s) \dg8.\ {t{s),q{s)) for a.a. s G (0,5), 

holds. A pair (t,q) : [so,Si] — > [0, T] x X satisfying 01.61) (with [0,5*] replaced by [so,Si]) 
is called parametrized metric solution of the rate-independent system (X, d, £). 

Indeed, the very focus of this paper is on getting insight into the properties of para- 
metrized metric solutions and comparing them with the other solution notions for rate- 
independent evolutions. That is why, in order to avoid technicalities and to highlight, 
rather, the features of our approach, throughout the next sections we shall work in a 
technically simpler setup, in which the state space X is a finite-dimensional manifold, en- 
dowed with a (Finsler) distance d associated with a 1-homogeneous dissipation functional 
D^ : TQ —>■ [0, cxd), and an energy £ G C^([0,T] x Q). The fully general metric framework 
is postponed to the forthcoming paper [MRS08J . see also Section [61 

The notion of parametrized metric solution [t, q) generalizes the outcome of the finite- 
dimensional vanishing viscosity analysis of |EfM06j and hence allows for the same mechan- 
ical interpretation, see Remark l3.5[ Namely, according to whether either of the derivatives 
f or Ig'l is null or strictly positive, one distinguishes in fll.6b|) three regimes: sticking, rate- 
independent evolution, and switching to a viscous regime (in correspondence to jumps of 
the system from one metastable state to another). In this metric setup as well, we show 
that the behavior of the system along a jump path is described by a generalized gradient 
flow. This can be seen more clearly when considering the non-parametrized solution q 
corresponding to the pair (t, q). The latter functions are called BV solutions of (Q, d, £) 
and are pointwise limits of the un-rescaled vanishing viscosity approximations g^, see Def- 
inition 14.31 and Section H] for an analysis of their properties. In particular, we shall show 
how to pass, by means of a suitable transformation, from a (truly jumping) BV solution 
g to a ("virtually" jumping) parametrized solution (t,q), and conversely. 

In Section [5], we compare the notion of BV solutions with other solutions concepts, 
namely with the energetic solutions of |MiT99| IMTL02J , and with the approximable and 
local solutions of |KMZ07i IToZ06[ |Cag08| (suitably rephrased in the metric setting, see 



Definitions 15. H 15. 5[ and 15.61) . In Section [^751 we review the notion of ^-minimal solutions 
of a rate-independent evolutionary system, proposed in |Vis01] using a global variational 
principle in terms of a suitably defined partial order relation between trajectories. First, 



we conclude that the notion of local solution is the most general concept, including en- 
ergetic and BV solutions, whereas BV solutions encompass approximable and ^-minimal 
solutions. Moreover, our notion of BV solutions has "more structure" , which makes it ro- 
bust with respect to data perturbations (cf. Remark l3.10p . whereas neither approximable 
nor ^-minimal solutions are upper-semicontinuous with respect to data perturbations. 

Further insight into the comparison between the various solution notions is provided 
by the examples presented in Section [TJ which are one or two-dimensional, such that the 
set of all solutions can be discussed easily. The one- dimensional case in fact relates to 
crack growth (under the assumption of a prescribed crack path), which was treated in 
[ToZOGt Cag08[ INeO07l IKMZ07] . The solution concepts developed there are also based on 



the vanishing viscosity method. In the latter case, the solution type in fact coincides with 
our notion of BV solution. We postpone the more difficult PDE applications to |MRS08j , 
where we are going to combine the notions of this paper with the methods of [RMSOSj 
to develop the present ideas in the infinite-dimensional or fully metric setting. Related 
ideas using the vanishing viscosity method for PDEs are found in |DD*07j . for a model 
for elastoplasticity problems with softening, and in |MiZ08j . for general parabolic PDEs 
with rate-independent dissipation terms. 

2 Setup and mechanical motivation 

We consider a manifold Q that contains the states of our system. The energy £ of the 
system depends on the time t G [0, T] and the state g G Q. Throughout the paper, we 
shall assume that £ G C^(Qt), where Qy = [0,T] x Q denotes the extended state space. 
The evolution of the system is governed by a balance between the potential restoring 
force — Dg£(t,g) and a frictional force /. The latter is given by a continuous dissipation 
potential D? : TQ -^ [0, oo), in the form / G dq'}l{q,q). We generally assume that 
Di(g, ■) : T^Q -^ [0, cxd) is convex and dq'Jl^q.q) C T*Q is the set-valued sub differential. 
Hence, the system is governed by the differential inclusion 

G dq-Riqit), eq{t)) + D/(i, ^W) C T^Q , t G (0, T) , (2.1) 

in which we have introduced a small parameter £ > to indicate that we are on a very 
slow time scale. 

Further, we suppose that 3^ = CRi+3^2, where 3^i : TQ -^ [0, oo) and 3^2 : TQ — > [0, oo) 
are such that for every g G Q 

CRi(g, ■) is convex and homogeneous of degree 1, 
^2{q, ■) is convex and homogeneous of degree 2. 

Note that 'Jlj{q,'yv) = ■y^^{q,v) implies d^j{q,'yv) = '-y^~^d^j{q,v) for all 7 > and 
(g, v) G TQ. Hence, 02. ip takes the form 

G dqJii{q{t), m) + sdqMQit), m) + D/(^, qit)) , t G (0, T). (2.2) 

We call 3^1 the potential of rate-independent friction and 'JI2 the potential of viscous 
friction. 



Remark 2.1 A prototype of the mechanical situation we aim to model arises in con- 
nection with a system of k particles {k > 1) moving in W^, hence with state-space 
Q = {q = (gi, . . . ,qk) : gj G W^} = ISJ"^ . We impose rate-independent friction ^i and 
viscous friction JI2 via 



^i(g,g) = Xl^^^i)l^jl ^^^ ^2(g,g) = XI^^I^jI' ^^^ (^'^) ^ 



fed X M'^'^' . 



where \qj\ is the Euclidean norm of the j-th particle velocity and /x, z/ : M'^ ^ [0, C)o) 
are given continuous functions. For k = 1 the potentials CR^ are related by CR2(g, g) = 
2^4fT ^i(q'; q), while for k = 2 their interplay is more complex. 

Our aim is to understand the limiting behavior of the solutions to (12. 2p for e ^ 0. In 
fact, we expect that for £ — > the rate- independent friction dominates, but the solution 
g^ : [0,T] — » Q may develop sharp transition layers, with q of order 1/e. In the limit we 
obtain a jump, but in order to characterize the jump path the viscous potential is crucial. 

The key idea is to study the trajectories T^ = { (t, g'^(t)) | t G [0, T] } in the extended 
state space Qt- The point is that the hmit of trajectories % may no longer be the graph 
of a function. To study the limits via differential inclusions, we may reparametrize the 
trajectories 7^ in the form 

7e = {iUs)Ms))\se[o,S,]}, 

where ts is supposed to be nondecreasing and absolutely continuous. 

For passing to the limit it is now helpful to select a family of parametrizations via 

m^ G Lj'q^((0, 00)) converging to m in Ll^^{{0 , 00)) , with ■m{s),m^{s) > for a.a. s G 
(0, 00), and to assume 

?,{s) + V^^S(s)Ms)) = me{s) for a.a. s G (0, ^,) . (2.3) 

Note that this can always be achieved. In particular, when 

Q = R^ 3l2(g,g) = -|gp V(g,g) gM'^x M"^ and m = m^ = 1 , (2.4) 

relation (12. 3p leads to the arclength parametrization of T^, which was considered in 
[EfMOBj . The total length is S"^ := T + J^ sj2'R2{qe{t),q'M) dt. Since ^^ ^ ^ up to 
a subsequence (thanks to standard energy estimates), it is not restrictive to assume that 
Se is independent of e by the simple linear rescaling fhe^s) = ms{sSs/S). 
By the chain rule and the j-homogeneity we have 

d^Ji,iqs{t),qemt=ns) = (p^) d^^RjiUs)^^)) fora.a. s G (0, 5,) . 
Now, using (12. 3|) and defining 



^{q, v) = gi^2iq, v)) with g{r) 



losll^Ti;)-^ forrG[0,i), 
00 otherwise. 



easy computations (cf. [EfMOGl Thm. 3.1] in the particular case of (12.41) ) show that (12.21) 
is equivalent to 

e d^Mq^q'^) + ed^dlU,^) +D,0„q,) a.e. in {0, S,) . (2.5) 

In this formulation we may pass to the limit, and we expect to obtain the limit problem 



Oed^Jl{q,l) + D,0,q), 



?+^2%{q,q 

where 

^{q,v): 



^^ — m. 



a.e. in (0,5), (2.6) 



Ml,v) ioi5l2{q,v)<l, 



oo for Ji2{q,v) > | . 

Formulation (12.61) is in fact a generalization of the one in |EfM06] . where rigorous con- 
vergence proofs of problem (12. 5p to (12. 6p are derived (in the case Q is finite-dimensional). 
Although 01 is no longer 1-homogeneous, the limit problem is still rate-independent: upon 
adjusting the free function m, one sees immediately that system (12. 6p is invariant under 
time reparametrizations. 

In the present work we concentrate on the case ^2{q,v) = |CRi(g, t>)^, since it is 
this case which can be generalized to abstract metric spaces and hence to the infinite- 
dimensional setting, see [MRSOSj . By introducing the dual norm of a co- vector w G T*Q 

Jii^,{q,w):=snp{{w,v):veTgQ, Jl^{q,v)<l}, (2.7) 

the operators d5li{q, ■) and d^2{(l, ■) can be characterized by 

w edJii{q,v), Vy^O ^^ Jli^,X(l,w) = l, {w,v) =Jii{q,v)>0, (2.8a) 

wed3ii{q,0) ^^ 3^i,,(g,w) < 1, (2.8b) 

wed%{q,v) ^^ Oli^,{q,w) = Jii{q,v) = {w,v), (2.8c) 

and they satisfy d^2{(l,'v) = 3^i(g, f ) 93li(g, f ) and 

'Jii^^{q,w)'Jli{q,v) = {w,v) <^^ w E \d'Jii{q,v) for some A > 0. (2.9) 



Proposition 2.2 In the case ^2{q,v) = ^Jii{q,vy, a pair (t,q) E AC([0, 5]; [0,T] x Q) 
fulfils (12. 6p (for some m G L^{Q,S) with m{s) > a.e. in (0, 5) j if and only if there 
exists a function A : (0, S) — > (0, oo) such that 

0GA93^i(g,g') + D,£(t;g), 

?>0, A > 1, (A-1)? = 0, ) a.e. zn (0,5). (2.10) 

? + aii(g,g') >0 



Proof: First, we note that, using the l-homogeneity of "Jii, the second of (12. 6p may be 
rewritten as ^ 

— + 0^1 (g,^ 1=1 a.e. in(0,5). (2.11) 

m \ mj 

Now, it is not difficult to see that 



a,5(f. i) = I ^''^'ffi a , 'f ^'(«- 1) ^ l°- 1) (-/ > °) ■ (2.12) 

" '^'m' \ |l,oo).S,S,(9,a if*,(f,J;) = l («i' = 0), 

where the equivalences in parentheses follow from (12. lip . Combining (I2.12p with the first 
of (12. 6p and using that d^i is 0-homogeneous, we deduce (I2.10p . 

Conversely, starting from (I2.10p . we put m{s) := t'(s)+CRi(g(s), q'{s)) for a.a. s G (0, S) 
and note that, by the third of (12.100 . m > a.e. in (0, S) and m G L^{0, S), since Oii is 
continuous. Using (12.110 and arguing as in the above lines, one sees that, if the pair (t , A) 
satisfies the second of (I2.10p . then Xdq^i{q,q') = dg^{q, ^) , which allows us to deduce 
the differential inclusion in (12. 6p from the one in (12.100 . ■ 



3 Analysis with metric space techniques 

3.1 Problem reformulation in a metric setting 

First of all, we complement the setup of Section [2] by specifying our assumptions on the 
rate-independent system (Q, 3?i, £), where Q is the ambient space, Jii the dissipation func- 
tional, and £ the energy functional. The more general setup will be studied in |MRS08] . 
Namely, we require that 

Q is a finite-dimensional and smooth manifold, (3.Q) 

and the energy functional satisfies 

£gC^(Qt). (3.E) 

The dissipation functional 3^i : TQ — ^ [0, oo) is a complete Finsler structure on Q (see e.g. 
[BCSOOl Ch.I.l]), namely 

5?i is continuous on TQ and Vg G Q : 3^i(g, ■) is a norm on T^Q , (3.R) 

called Minkowski norm in the Finsler setting. Then, CRi induces the (Finsler) distance 
rf: Qx Q^ [0,oo): 

d{qo,qi) ■■=mm^ / 5ii{q{s),q{s))ds:qeA{qo,qi)j, 

where for all qo,qi & Q we set 

MQo, qi) = {ye AC([0, 1], Q) 1 1/(0) = go, 2/(1) = Qi }• (3.1) 

8 



Hence, (Q, d) is a metric space, which we assume to be complete. Like in the previous 
section, we let 3^2 = ^^i- For a curve q G AC([0, T]; Q), the Finsler length of its velocity 
q'{t) is given by 

\q\{t) := 5ii{q{t),q{t)), well-defined for a.a. t G (0,T), (3.2) 

and satisfies 



d{q{s),q{t))< \q\{r)dr for all < s < t < T . 
Using Jli, we define the associated local slope |(9g£| : Qr — ^ [0, oo] via 



(3.3) 



|a,£|(t,g):= sup iM(^l^ = 3?i,,(g,D,8(t,g)), (3.4) 

which is the conjugate norm with respect to the Minkowski norm 3^i(g, ■) of the differential 
of the energy in the cotangent space T*Q. 

Using the smoothness of £ we have that for every curve (t,q) G AC([so,Si]; Q^) the 
map s ^— s> £(t(s), q{s)) is absolutely continuous and the chain rule for £ gives 

^£(t(s), q{s)) = dMt{s), q{s)) t'{s) + (D,£(t(.), q{s)), q\s)) (3.5) 

for a.a. s G (sq, Si). On the other hand, formulae (13. 2p and (13.41) yield 

(D,£(t,g),g'))>-|g1 \d,^\{t,q) . (3.6) 

Therefore, every (t, q) G AC([so, Si]; Qt) fulfills the chain rule inequality 

^^e{t{s),q{s)) - dtms),q{s))t'{s) > -|9,£| (t(s), g(s)) |g'|(s) a.e. in (sq, Si). (3.7) 

The metric formulation of doubly nonlinear equations. We now see how notions 
(13. 2p and (13. 4p so far introduced come into play in the reformulation of a class of doubly 
nonlinear evolution equations in the metric setting {Q,d). 

Let ip : [0, oo) -^ [0, oo] be a lower semicontinuous, nondecreasing, and convex function 
and ip* : [0, oo) -^ [0, oo] its conjugate function (Legendre-Fenchel transform), namely 

7/>*(O = sup{z/^-^(z/)|z/>0}. 

Following |AGS05j (see also |RMS08] ). a function q G AC([0,T]; Q) is called a solution of 
the '^/'-gradient system associated with (Q, (i, £) if 

^£(t, q{t)) < dMt, q{t)) - ^iWm) - r{\d,8-\ (t, q{t)) ) a.e. in (0, T) . (3.8) 

It has been proved in |RMS08[ Prop. 8.2] that q fulfills (13. 8 p if and only if it solves the 
doubly nonlinear equation (also called quasi- variational evolutionary inequality) 

G dg^{q{t), q{t)) + D,£(t, q{t)) a.e. in (0, T) , where ^(g, q) := V^(3?i(g, q)) . (3.9) 



Under assumptions (|3.QD, ( 13. Rp . and ( 13.Ep . the existence of absolutely continuous solu- 



tions to the Cauchy problem for fl3.8p follows from |RMS08t Thm. 3.5]. We stress that 
the simple, but central duality inequality '?/'(//) + V'*(0 ^ ^^ fo^' all i^, ^ & [0, oo), together 
with the chain rule inequality (13.71) . enforces equality in (13. 8 p (ultimately in (13. 7p as well). 
In the rate-independent setting, the natural choice is 

^o{iy) = V giving ^^(0 = /[o,i](0 > 

where /[o,i](') denotes the indicator function of [0,1], i.e. /[o,i](0 = if ^ G [0,1], and 
-^[0,1] (0 — ^^ otherwise. However, simple one- dimensional (not strictly convex) examples 
show that we cannot expect existence of absolutely continuous solutions in this case, cf. 
Example 17.11 Hence, we proceed as in Section [2] and consider limits of viscous regulariza- 
tions after suitable reparametrizations. 

Before doing so, note that (13. 8p is equivalent to the parametrized version on some 
interval (sq, Si), given by 



££(f(.),g(.))-9,£(f(.),g(.))?(.) 



(3.10) 



where g(s) = ^(hs)) and t'(s) > a.e. in (so,Si). In the rate-independent case, the 
right-hand side does not depend on t (s), because '^^iy) = v implies i\)Q{av) = aipo{i') 
and ipoiO = (^i^oiO ^^^ all a > 0. 

3.2 Rate-independent limit of viscous metric flows 

We now consider the case of small viscosity added to the rate-independent dissipation, 
namely 

^^(z/) = z/+£z/2 Vz/e[0,oo). (3.11) 

We obtain ^*(0 = ^((^-1)+)^- Thus, (IXTUD takes the form 

E{t{s),q{s))-dMt{s),q{s))^is)<-M,{t'{s),\q'\{s),\d,E\{t{s),q{s))) (3.12) 



2e 
dl 

for a. a. s G (sq, Si), with 



M,(«,^,e):=«^.(9+«V^:(0 = ^+^^' + ^((e-in' (3.13) 

for all (a, u, ^) G (0, cxd) x [0, oo)^. Clearly, for fixed a > the limit as e \ of M^{a, u, ^) 
gives ipo{i') + 'ipoiO- However, our purpose is to blow-up the time parametrization when- 
ever jumps occur. Indeed, the finite-dimensional case (see |EfM06j ) suggests that jumps 
in the rate- independent evolution will occur at fixed rescaled time (i.e., when t' = 0). 
Hence, we also have to consider the case a ^ as e —* 0. For this, note that when ^ > 1 
Mir^-jUyC,) assumes its minimum on [0, oo) for al = eu/{^ — 1)+, corresponding to the 
value M^{a.l, u,^) = z/ + z/(^— 1)+. In any case we have 

M,(a,z/,e) >Minf(z/,0 :=z/ + z/(^-l)+ for all (a,z/,0 e (0, oo) x [0,oo)^ (3.14) 
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Thus, we define Mq : [0, oo)'^ — > [0, oo] via 

M„(a...O := I *'■"''"■«> = ''^"''«-'>' '"" = "■ (3.15) 

1 Af,„p((/,f) = >'+/MK) fora>0, 

and obtain the following result. 

Lemma 3.1 Define M^ : [0, oo)^ -^ [0, oo] wa fIXTSD . M^(0,0,O = for all ^, and 
M^{0, z/, ^) = oo for all ^ and u > 0. Then, we have the following results: 

(A) Mfr T -converges to Mq, viz. 



ae,Ue,^e) -^ ia,^,0 =^ Mo{a,u,^) < liniinf Me(a£, i/e,^^ 

£\0 



(3.16a) 



T-liminf estimate: 

T-limsup estimate: 

w. c\^(( c\\ \ {ae.Pe.ie) ^ {a,p,i) and (3.16b) 

V(a,i/,0 3((a„i/„^,)),>o: \ 

[ Mo(a, u, 4) > hmsup^^o ^e(«£, ^e, 4e) • 

(B) If {ae.Ve) -^ (5?, z/) in L-'^((so, Si)) and liminfe^o^e('5) > i{s) a.e. in (so,Si), then 
Mo(a(s),i?(s),^(s))ds<liminf / M,{aeis),iy,{s),Us))ds . 

Proof: Estimate (I3.16ap is trivial for a > 0, as we have pointwise convergence then. If 
a = 0, we employ (13. lip and use that Mjnf is continuous. 

To obtain (]3.16bp in the case a > we simply take (a^, i^e,^^) = {a,^,^) and the 
result follows from pointwise convergence. If a = 0, we let (a^, u^, ^s) = («*, ^, ^^^ the 
desired result follows. Thus, (A) is proved. 

To show the estimate in Part (B), let us introduce the function M : [0, oo)^ -^ [0, oo] 

M(a,z/;e,£):=Me(a,z/,0; 

by the previous point (A) and the fact that M is lower semi continuous when e > 0, it 
is immediate to check that M is lower semicontinuous in [0, oo)^. Moreover, M (-, ■; ^,5) 
is convex in [0, 00)^ for all C,,b: this property can be directly checked starting from the 
definition of M or by observing that Ms{-, ■, ^) is convex when e > thanks to (I3.13P and 
the convexity of the map (z/, a) 1— > u'^/a. 

Assuming initially that ^^ — >■ ^ in L^{si, S2) and considering an arbitrary infinitesimal 
subsequence e„ — >• 0, we can then apply loffe's Theorem (see [lof 77] ) to the sequence of 
maps s ^ («£„(«), z/e„(s),^e„(s),5„), obtaining 

f ' Mia{s),i}{s),as),0)ds<\immi /" ' M(«,„(s), z/,„(s),e.„(s),£„)ds . 

/ n— »oo / 

■J so J so 
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In the general case, we consider an arbitrary k > and we replace ^^^ with the sequence 
^^,n{s) := min(^£„(s),^(s),K), converging to ^^(s) := min(^(s),K) in L^(si,S2). Since M 
is nondecreasing with respect to ^, we argue as above and obtain 



n - 
M(a(s), i/(s), ^«(s), 0) ds <liminf / M(a£„(s), i/£„(s),^«,„(s),5„)ds 

< lim inf / M{a,^ (s) , u,^ (s) , ^,^ (s) , e„) ds . 



Sl 



Passing to the limit as k ^ oo and applying Fatou's Lemma we obtain the desired 
inequality. ■ 

In fact, the specific form of Mq is not needed in the sequel. Hence, we will consider 
general functions M : [0, oo)^ -^ [0, oo] with the following properties (which are obviously 
satisfied by Mq): 

M : [0,oo)^ ^ [0,oo] isl.s.c, (3.17a) 

M(7a,7z/,0 =7M(a,i/,0 for all a, i/,^,7 G [0, cx)) , (3.17b) 

M{a,u,^)>u^ for alla,i/,^G [0,oo), (3.17c) 

M(«,z/,e)=z/e ^^ («,z/,OeH, (3.17d) 

where the set S := S^^^'^^ u H'^^'p U 3^™^ consists of the disjoint fiat pieces (see Figure IXTI) 

S^tick _ I (^^ 0^ ^) ^ [0, oo)3 I a > 0, e G [0, 1) }, 

S^iip := {(a,i/,l) G [0,oo)3|a >0,z/>0}, and (3.18) 

^jump ._ I ^Q^ j^^ ^^ G [0, oo)=^ I z/ > 0, ^ > 1 }. 

For instance, the function 

M{a, u,^) = iy+ (^-l)+(a+z/) = max(^, l)z/ + (^-l)+a (3.19) 

fits in this framework. It is not difficult to check that, if M is nondecresasing with respect 
to ^, then M < Mq. 

The notion of parametrized metric solutions. The following notion of parametri- 
zed metric solution of the rate-independent system (Q, d, £) is proposed in a general form, 
replacing the function Mq obtained in the vanishing viscosity limit with a generic function 
M satisfying (I3.17p . The proposed notion is fitted to the metric framework and does not 
need a differentiable structure. However, it strongly relies on the fact that the small vis- 
cous friction £:|g'|(t)^ is given in terms of the same metric velocity as the rate-independent 
friction, see also the assumption 3^2 = |3^i- We refer to [EfMOGl IMiZOS] for settings 
avoiding this assumption. 



Definition 3.2 (Parametrized metric solution) Let {Q,d,8,) satisfy p.QP, (13. Rp . and 



(13. Ep and let M fulfill (I3.17p . An absolutely continuous curve {t, q) : (sq, Si) -^ Qt is called 
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Figure 3.1: The set E = S"*'^'^ U E^^^'p U 3""^? 



a parametrized metric solution of (Q, d, £), if 

t : (so,Si) -^ [0,T] is nondecreasing, 



t (s) + \q\{s) > for a. a. s E (sq, Si), 

^^E{t{s),q{s))-dtms),q{s))t'{s) 
< -M{?{s), \q'\{s), |9,e| {t{s), q{s)) ) /or a.a. s G (sq, s^ 



(3.20a) 
(3.20b) 

(3.20c) 



If {t,q) satisfies only (I3.20ap and (I3.20cl) . it is called a degenerate parametrized metric 
solution. Ift : (so,Si) — >• [0, T] zs also surjective, i.e. t(so) = and t{si) = T, then (t,q) 
is called a surjective parametrized metric solution. 

This solution concept has the concatenation property as well as the restriction property. 
The former means that if (t,g) : (sq, Si) — >• Qt and (t,g) : (si, S2) — ^ Qt are parametrized 
metric solutions with (t(s7), g(s^)) = (^1,^1) = (t{sf),q{sf)), then their concatenation 
(t, g) : (so, S2) ^ Qt is a solution as well. We point out that, thanks to fl3.17bp . the 
notion of parametrized metric solution is rate- independent, i.e., invariant under time 
reparametrizations by absolutely continuous functions with strictly positive derivative a.e. 
(by nondecreasing absolutely continuous functions in the case of degenerate solutions). 
Moreover, the notion is independent of the particular choice of M , as long as M satisfies 

(EIZD- 

Remark 3.3 (Nondegeneracy and arclength reparametrization) Any degenerate 
parametrized metric solution admits a nondegenerate reparametrization (t, q) : [0, S] — >■ 
Qt, thus satisfying also (]3.20bp . This means that t(s) = t(cr(s)),g(s) = q{a{s)) for some 
absolutely continuous, nondecreasing and surjetive map a : [sq, Si] -^ [0, S]. In particular, 
we can choose a so that i' + \q'\ = 1 a.e. in (0, S) by defining 

ais):= r {?(s) + \q'\{s))ds = t{s) ^t{so) + [' W^ds, S := a{s^) 



so 



so 
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(cf. also Lemma [4. ip . In fact, for every interval [rojTi] C [sq, Si] we then have 

or(ro) = o-(ri) -^ f(ro) = f(r) = ?(ri), g(ro) = g(r) = g(ri) for all r G [ro, ri]. 

We can then define {i{a),q{a)) := (t(s),g(s)), whenever a = a{s). For ao = cr(ro) < ai = 
cr(ri) we obtain 

t(cTi) -t(cro) + c?(g(o-i),g((To)) < / {?{s) + \q\{s))ds = ai - (Xo , 

Jro 

giving t' + |g'| < 1 a.e. in [0, S]. The nondegeneracy condition holds with t' + |g'| = 1 a.e. 
in S, which follows via the change of variable formula: 

^> [\i' + \q'\)da= r {i'iais)) + \q'Ms)))a'is)ds 

Jo Jso 

'so 

Parametrized metric solutions admit various different but equivalent metric charac- 
terizations (where we avoid to explicitly use the differential Dg of the energy). 



{?{s) + \q'\{s))ds = a{s,) = S. 



Proposition 3.4 Under the same assumptions of the previous Definition l^/.j^ o't- abso- 
lutely continuous curve (t,q) : (scSi) -^ Qt satisfying fl3.20ap and fl3.20bp is a parame- 
trized metric solution of (Q, d, £) if and only if it satisfies one of the following conditions 
(equivalent to fl3.20cp ).• 

A) For all sq < ctq < cxi < Si we have 



£(f(ai),g(ai)) - £(f(ao), g^) - T dte(t{s),q{s))?{s) ds 

J (TO 



(3.21) 



<- / M{?is),\q'\is),\d,E\{tis),qis)))ds. 



O-Q 



B) Eqn. (13.211) holds just for ctq = Sq and ai = Si, i.e. 



Sl 



£(t'(si),g(si))-£(t(so),g(so))- / dtE{t{s),q{s))?{s) ds 

J so 

< - r M{?{.s), Ig'l(s), |9,£| {tis),qis)) ) ds. 



(3.22) 



IVI "^ 
J so 

C) For a. a. s G (sq, Si) we have 

^^E{t{s),q{s))-dt&{t{s),q{s))?{s) = -\q'\{s) |9,£| {t{s),q{s)) , (3.23) 

and one of the following (equivalent) properties 

M{?{s),\q'\{s),\d,8\ {t{s),q{s))) = |g'|(s) |a,£| {t{s),q{s)) , (3.24a) 

{?{s),\q'\{s),\d,^\{t{s),q{s))) G S, (3.24b) 

?(s)>0 =^ \dgE\{t{s),qis)) <1, 
|g'|(,)>0 =^ |9,£|(f(s),g(.))>l. 
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(3.24c) 



In particular, by (13.231) and f l3.24ap . the following identity holds a.e. in (so,si); 

^£(f(s),g(s)) - dtms),q{s))?{s) = -M{?{s), \q'\{s), |9,£| {t{s),q{s))). (3.25) 

Proof: A) is just the integral formulation of (I3.20cp . 

We note that the chain rule inequality (13.71) . combined with (I3.20cl) and (I3.17cp . implies 
(I3:23|l and ^2M . By condition (l317dl) . (jSJUD is equivalent to ^ME^ . 

Since the set S can be easily characterized via 

{a,u,^)eE ^^ ("(a >0 ^^ < 1) and (z/ > 0^^> 1) 

we ultimately find that fl3.24bl) can be replaced by the simple relations (I3.24cl) . 

Having obtained the equivalence between (I3.20cl) and C), we can now show that B) 
is sufficient to characterize parametrized metric solutions (the necessity is trivial): again 
applying the chain rule (I3.5p - (I3.7I) . we get 

Tit 



M{?{s)^q!\{s)^^,^\ {t{s),q{s))) - (-D,£(t(s), g(s)), g'(s)) J ds < 

so that fl3Tra and ([M]) yield 

M(?(.),|g'|(s),|9,£| (t{s),q{s))) = (-D,£(f(.), g(.)), g'(.)) = |g'|(.) |9,£| (f(.),g(.)) 
for a. a. s G (so,Si). We thus get (13.231) and (I3.24ap . ■ 

Remark 3.5 (Mechanical interpretation) The evolution described by relations (I3.24cp 
bears the following mechanical interpretation, cf. |EfM06] . Indeed, with (a, z/, ^) = 
(?, Ig'l, |5g£| {t,q)) we can use the decomposition S = E^^^'^^ u S^^''^ U SJ"™^: 

• (? > 0, Ig'l = 0) leads to sticking {{a, z/,0 e E""^'""^), 

• (f > 0, Ig'l > 0) leads to rate-independent evolution ((a, i/, ^) G 2*^''^), 

• when [t' = 0, Ig'l > 0), the system has switched to a viscous regime, which is 
seen as a jump in the (slow) external time scale (the time function t is frozen and 

(a,z.,0e3"-p). 

Remark 3.6 Properties (I3.17cp and fl3.17dp seem to be related to the notion of bipotential 
(cf. e.g. |BdV08j ). which was proposed for studying non-associated constitutive laws in 
mechanics by convex analysis tools. We recall that, given two (topological, locally convex) 
spaces in duality Z and Z', a function 6 : Z x Z' — > (— oo, oo] is called a bipotential if 
it is convex, lower semicontinuous with respect to both arguments, and fulfills for all 
(z/, ^)eZ X Z' 

b{u,o>{^,^) and (eG9,6(-,0M ^ i^ed^b{u,-m ^ K^,o = (e,^) 

Indeed, for every a > the functions Me(a, -, ■) given by (13.130 and Mo(a, ■, ■) by (13.150 
are bipotentials on [0, oo) x [0, oo). 
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The next result ensures that the abstract metric evolution formulation developed here 
reduces to the one stated in Proposition I2.2[ 



Proposition 3.7 Let {Q,d,8.) satisfy p.QP, (13. Rp . and (13. ED . Then, a curve (t,q) G 



AC([so, si]; Qt) is a parametrized metric solution of the rate-independent system (Q, d, £) 
if and only if there exists A : (sq, Si) —>■ [1, cxd) such that (12.101) holds. 

Proof: By (Q, condition (13.20bD in Definition O coincides with the third of (12101) . 
Now, let us first suppose that (13.20cp holds, and set A(s) := max{|c)g£| (t(s),g(s)) , 1} . 
We shall prove that the triple (f,g,A) fulfills (l^TTUj) on (so,Si). Indeed, (E^SD, (IS3D, and 
(El]) yield 

- D,£(f(s), g(s)) G |9,£| (t(s), g(s)) 9D?i(g(s), g'(s)) for a.a. s G (sq, si) . (3.26) 

Now, let us fix s G (sq, Si) at which (13.26P holds: if \q'\{s) > 0, taking into account 
the second of (13.24cp we find that A(s) = \dqE,\ (t{s),q{s)) and that, by (13.26p . the triple 
(r, |g'|. A) satisfies the first of (12.10p at s = s. On the other hand, if \q'\{s) = 0, necessarily 
?(s) > by (13.20bp and the first of (KMch gives that \dgE\ (f(s),g(s)) < 1. In this case, 
A(s) = 1 and (13.260 implies 

-D,£(f(s),g(s))Gaaii(g(s),0), 

hence we again conclude that (t'{s), |g'|(s), A(s)) fulfills (12.10p ^. 

Conversely, from the first of (12.10p we read that for a.a. s G (sq, Si) 

\d,E\ {t{s),q{s)) Wl = {-DgE{t{s),q{s)),q'{s)), (3.27) 

\d,E,\{t{s),q{s))<X{s), (3.28) 

\dqE\ (f(s),g(s)) = A(s) if \q\ > 0. (3.29) 

Then (IX^UcD follows, if we check (K^^ . Indeed, if ? > 0, then the second of d^TTU]) yields 
A = 1. Therefore the first condition of (13.24cp follows from (13.280 . If \q'\ > 0, combining 
( 13.29P and the constraint A > 1 of (12.100 . we also get the second of (13.24cp . ■ 

Convergence of the viscous approximation. The main result of this section states 
that limits (t, q) of parametrized solutions (t^, qe) of the viscous system (13.100 . with ip = 
ipe, are actually parametrized metric solutions of the rate-independent system (Q, (i, £). 
By the standard energy estimates and an elementary rescaling, it is not restrictive to 
assume that the domain of (t^, q^) is a fixed interval [sq, Si], independent of e. 



Theorem 3.8 (Vanishing viscosity limit) Let(Q,d,E,) satisfy (|3.QP, ( 13. Rp . and ( 13. EO 



For every e > /et g^ G AC([0,T]; Q) be a solution to (13. 8p for ip = Tp^. Choose non- 
decreasing surjective parametrizations is G AC([so,si]; [to,e,T]), and let qe{s) = ge(te(s)). 
Suppose that there exists go ^ Q, o^nd m G L-'^((0, S)) such that 

to,6=%iso) ^0, qei-so) = qeito,e) -^ Qo ase\0, (3.30) 

me := t^ + \qj ^^ m inL^(so,Si) ase\0. (3.31) 
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Then, there exist a subsequence ((t^^, g^^))fcgN with Sk \ and (t,q) G AC([so, si]; Qt) 
such that Ct{so),q{so)) = (0,go); '^'^t?; as k -^ oo, 

{%^,qj^{t,q) znC°([so,Si];QT)), (3.32) 

{t's,,WsJ)-it',W\) ^nL\[so,s,]■,R'), f \q'J{t)dt ^ T \q'\{s)ds. (3.33) 

The limit {t,q) is a degenerate parametrized metric solution of {Q,d,C,) (i-e. it satisfies 
(I3.20al) and (I3.20cl) ). and it is nondegenerate (recall fl3.20bl) ) if m{s) > a.e. in {sq,si). 

Proof: Eqns. (EJD and (l33B yield 

diqs{ro)Mri))< /V.l(5)ds< f'm,{s)ds. (3.34) 

In particular, choosing Tq = Sq and using fl3.30j) - fl3.31j) we find C > such that 

d{qo, Qeit)) < C for all t e [sq, Si] and all e > . (3.35) 

Moreover, it follows from fl3.3ip that the sequences {t^} and {|g^|} are bounded and 
uniformly integrable in L^(so, Si). Hence, on the one hand, the Ascoli-Arzela compactness 
theorem and its version for metric spaces |AGS05l Prop. 3.3.1] yield that there exists 
an absolutely continuous curve (t,g) : [sc-Si] — >• Qr such that, up to a subsequence, 
convergences fl3.32p hold. On the other hand, by the Dunford- Pettis criterion (see, e.g., 
[DuS88l Cor. IV.8.11]), there exists rj E L^{sq, si) such that, up to the extraction of a (not 
relabeled) subsequence, 

t^ ^ t , 15^1 ^^ ^7 inL^(so,si) as e \ 0. (3.36) 

Passing to the limit in (I3.34p we easily get 

\q\{s) < ri{s) < m{s) for a. a. s G (sq, Si) . (3.37) 



Now, the smoothness of £ and assumption ( |3.Q[ ) yield that £, (9t£, and \dgE\ are 



continuous with respect to both arguments, thus we readily infer from (I3.32p that 

06, Qe) -^ £(?,?), 1^/1 {te,qe) ^ \dg^\ {t,q) and dt8,(te,qe) -^ dtE,(t,q) 

[o.oo) 
uniformly in [sq, Si] for e \ . 

To proceed further, we integrate (I3.12p over [sq, Si] and obtain 

EiUso),Uso)) - 0e{si),qe{si)) + J 9i£(t;(s), g,(s)) ?,(s) ds 

° (3.39) 

> r M,(?,(.), %\{s), \dgE\ {Us), Us)) ) ds . 

J so 
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On the left-hand side we can pass to the hmit e ^ using (13.361) and f l3.38p . whereas for 
the right-hand side we use Part (B) of Lemma I3.lt 



mso),q{so))-msi),q{si))+ / dMtir) , q{r)) t\r) dr 



lim ( £(t,(so), g;(so)) - £(te(si), Usi)) + I dMt,ir),qe{r)) t[{r) dr 



Sl 



so 



> hniinf f ' M,(?,(r), \q'^\{r), \d,E,\ {Ur),Ur)) ) dr 

> r Mo{?{r),r]{r),\d,^\ (f(r), g(r)) ) dr 



J so 

> r Mo(?(r), |g'|(r), \d,8,\ (?(r), g(r)) ) dr 

J so 



where we have used (13.371) and the monotonicity of Mo(a, ■,^) for the last estimate. 

We see that (t, q) fulfills (13.221) a.e. in (sq, si) , and it is therefore a (possibly) degenerate 
parametrized metric solution of (Q,(i, £). Moreover, comparing the last inequalities with 
the integrated form of (I3.25p . we get 



Mo(?(r),r7(r),|9,£|(f(r),g(r))) = Mo(?(r), |g'|(r), |9,£| (f(r), g(r)) ) < 



oo 



for a. a. r G (sq, Si). Since Mo(a, -j^) is strictly monotone in its domain of finiteness, we 
get ri{r) = \q'\{r) for a. a. r G (so,Si), thus obtaining (I3.33p . Using the first convergence 
in (I3.36P we also find v + |g'| = m and the last assertion follows. ■ 

Remark 3.9 (Preservation of arclength parametrizations) If (t^, q^) are arclength 
parametrization (i.e. rri^ = 1), then their limit {t,q) still satisfies the arclength property 
? + Ig'l = 1, thanks to fl333D . This generalizes |EfM06l Cor. 3.6]. 



Remark 3.10 Mimicking the argument of the proof of Proposition 13.81 under the same 
assumptions it is also possible to prove a result of stability with respect to initial data 
for parametrized metric solutions. Namely, let {(t„,g5^)} be a sequence of parametrized 
metric solutions on a time interval [so,Si], such that (t„(so), gn(so)) = (^0)?o) ^^^ every 
n E N, with (^qj^q) -^ (^0, Q'o) and m^ := t'„ + |g|j| ^^ rri in L^{so,Si). Then, there 
exists a parametrized metric solution (too, 9oo), starting from (to, qo), such that, up to the 
extraction of a subsequence, (t„,g„) -^ (too,§oo) uniformly in [so,Si], with (t'^, |g^|) ^^ 
(t^Q, \q'ao\) in L^(so,Si). In fact, in [MRSOSj we shall prove the above result, as well as 
the vanishing viscosity analysis of Theorem 13.81 in the more general setting detailed in 
Section [61 

4 BV solutions 

Before introducing the notion of BV solution to the rate-independent system driven by 
£, we recall some definitions and properties of BV functions on [0, T] with values in the 
space (Q, d) introduced in the previous section. Note that, however, the following notions 
are indeed independent of the Finsler setting ( |3.QP - (]3.RP and can be given for a general 
complete metric space. 
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Preliminaries on BV functions. Given a function q : [0,T] ^ Q and an interval 
/ C [0, T], we define its variation on / by 

n 

Var(g, J) := sup ^ d{q{Tj^i), q{Tj)), (4.1) 

i=i 

where sup is taken over all n G N and all partitions tq < ri ■ • • < r„_i < r„ with tq, r„ G /. 
We set 

BV([0,T]; Q) = {g : [0,T] ^ Q I Var(g, [0,T]) < oo }, 

where we emphasize that functions are defined everywhere, as is common for rate-independent 
processes. For q G BV([0,T]; Q) and t G [0,T] the left and right limits exist: 

q{t~) := limg(t— /i) and qit~^) '■= \im.q{t+h), 

h\0 h\0 

where we put ^(0") = g(0) and q{T^) = q(T). In general, the three values qit~),qit), 
and q(t^) may differ. We define the continuity set Cg and the jump set Jg by 

Cg = { t G [0, T] I g(t-) = q{t) = q{t+) }, J^ = [0, T] \ C^. 

Indeed, our definition of "Var" is such that we have for all < r < s < t < T 

Var(g, [r, s]) = d{q{r), q{r+)) + Var(g, (r, s)) + d{q{s-), q{s)) , (4.2) 

and the additivity property 

Var(g, [r, t]) = Ya.i{q, [r, s\) + Var(g, [s, t]). (4.3) 

When calculating the variation of q over an interval J, one has to be careful with (possible) 
jumps at the boundary of /, if / contains boundary points. Now, for a function q G 
BV([0,T], Q) we introduce the nondecreasing function 

V, : [0,T] ^ [0,oo) , V,{t) := Var(g, [0,t]). 

The distributional derivative of Vq defines a nonnegative Radon measure /i^ such that 

/ig([s,t]) = Vqit) - Vq{s) Vt, s G Cg, (4.4) 

and, more generally (see |Fed69t 2.5.17]) 



/ at)dVq{t)= f CW/^,(dt) forallCGC°(0,T), 
Jo Jo 



(4.5) 



where J^ CdVg denotes the Riemann-Stieltjes integral. As usual, fig can be decomposed 
into a continuous (also called diffuse) part /x™ and a discrete part /x^, where for a Borel 

set y4 C [0, T] we have 

/xJ(A)=/.,(AnJg)= Yl diqit-),qit)) + diqit),qit+)), (4.6) 

teAnJq 
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in accordance with formula (14. 2p above. 

In the following technical lemma (whose proof is postponed to the end of this section), 
we will discuss the link between a BV map q : [0, T] ^ Q and its graph q(t) = (t, q(t)) in 
the extended state space Qr, endowed with the distance dQj.{{to,qo), iti,qi)) := |to— ^i| + 
d{qo,qi) ■ We denote by ^[a,b] the Lebesgue measure on the interval [a,b], whereas £ 
denotes a general one-dimensional Lebesgue measure. 

Lemma 4.1 Let q eBV{[0,T];Q) anc? q G BV([0,T]; Q^) with q{t) := (t, q{t)) . Set 

pit) ■= V,{t) = Ve,T{q, [0,t]), R := piT); a{t) := V^{t) = Var(q, [0,t]), S := a(T), 

with their right- continuous inverse 

f(r) := sup{ t e [0, T] | Vg{t) = p{t) < r }, f(s) := sup{ t G [0, T] | Vq(t) = a{t) <s}. 

Then, the following statements hold: 

A) a{t) = t+p{t), Jq = Jg, Cq = Cg, /iq = £ + /i„ /i™ = £ + /i™, /i^ = yU^ . 

B) There exist 1-Lipschitz maps q = (t, g) : [0, S] -^ Qy and p : [0, 5] ^ [0, R] such that 
q(t) = q(cr(t)), p(t) = p(o"(t)) /or all t G [0,T]. T/ie map t is uniquely determined, 
it is the right- continuous inverse of a, and it is infective on Cq := cr(Cq) = t^^(Cq). 
The maps q and p are uniquely determined on the set Cq and satisfy q{s) = q{t{s)) 
and r(p(s)) = t(s) for all s E Cq. 

C) t-^(£[o_5']) = Pq and r#(£j[o,/?]) = /^g; ^'^ ^^e sense that for all bounded Borel function 
C : [0, T] ^ M and Borel set A C [0, T] 



C(t(5))ds= / CWMdi), / C(r(r))dr= / CWp,(dt). (4.7) 

In particular, if A G [0,T],B G Cq G [0, S] are Borel sets, then 

p^{A) = L{r\A)), L{B) = /iq(f(i?)). (4.8) 

D) The Lebesgue densities of the measures £,Pg° ^ Pq° with respect to p^ are expressed 
by the formulae 



,co 



d-C ■:> ^f^q |-/| -7 



too, — = \q \ o a = p o a. (4.9) 



dp-° dp: 



The notion of BV solution. Let us first introduce a new family of 1-homogeneous 
dissipation functional So(t; ■, ■) : TQ -^ [0,oo), depending on the two parameters a G 
[0, oo) and t G [0,T], defined as 

§a{t', q, v) := max {|9g£| (t, q) , a} ^i{q, v) (4-10) 
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Notice that for all a > 0, t G [0, T], and g G Q the functional Sa(t; q, ■) is a norm on TgQ 
(possibly degenerate, when a = 0), thus satisfying condition (13.Rp . As in Section [XTl we 
can therefore consider the corresponding Finsler distances Sa(t; -, ■) : Q x Q — ;> [0, oo) via 



S a{t; qo, qi) := ini ^ / §a{t;y{s),y'{s))ds yeA{qo,qi] 



(4.11^ 



The functional 5*0 is called slope distance and Sa{t] ■, ■) also admits the equivalent formu- 
lation in terms of the metric velocity 

Sa{t;qo,qi) ■.= mi^ / meix{\dg8,\(t,y{s)) ,a] \y'\{s)ds y eA{qo,qi)j . (4.12) 

For a > the infimum in (14. lip and (I4.12p is attained. 

A straightforward consequence of the symmetry 3^i(g, — 1>) = 'Jli{q,v) (see (13. RP ) is 
that Sa(t;qo,qi) = Sa(t;qi,qo). Using the chain rule inequality (13. 7p we find 

|£(t,gi)-£(t,go)| < So{t;qo,qi) < S^{t;qo,qi) for all (t, go, gi) G [0,T]xQxQ. (4.13) 

The notion of BV solution to the rate-independent system (Q, d, £), which we are going 
to introduce, relies on a version of the chain rule for BV functions with values in a metric 
space. In order to state it, for a general q G BV([0, T]; Q) and < to ^ ^i ^ ^ we define 



So(g, [to,ti]) := 



\d,E\ (r, g(r)) /i^°(dr) + ^o(to; g(to), g(t+)) + ^o(ti; g(tr), Qih)) 



to 



+ Yl t^o(t; qin, qit))+Soit; g(t), g(t+))] . 
t&Jqn{to,ti) 



(4.14) 



Based on (I4.13p . we define a second functional T via 



r(g,[to,ti]):= / |a,£|(r,g(r))/i-(dr) 

Jto 

+ |£(to,g(to))-£(to,g(4))| + |£(ti,g(tr))-£(ti,g(ti))| 

+ Yl [^t,q{t)y8,{t,q{t+))\ + |£(t,g(r))-£(t,g(t))| 
tGJgn{to,ti) 



(4.15) 



Obviously, we have So(g, [s,t]) > T{q, [s,t]) > and both functionals So(g, ■) and r(g, ■) 
fulfill the additivity property (14. 3p . when considered as functions on intervals. 

Proposition 4.2 Under assumptions ( |3.QP , (13. Rp . and fl3.Ep . the following chain rule 
inequality holds for all q G BV([0, T], Q) and < to < ti < T.- 

£(ti, g(ti)) - £(to, g(to)) - / ' dt8,{t, q{t))dt > -T{q, [to, t,]) > -Eo(g, [to, t,]) . (4.16) 
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Proof: The function t \-^ E{t) := £(t, g(t)) is of bounded variation on [0, T] and its jump 
set is contained in Jg. We denote hy rj = ^E its distributional derivative (a bounded 
Radon measure on (0,T)) and by 7]'^° its diffuse part, defined as r]'^"{A) := r]{A fl C,) for 
all Borel sets Ac [0, T]. Thanks to (14.131) . we have fl4.16p if we show that 

r7^°>5,£(-,g(-))/:-|9,£|(-,g(-))^r (4-17) 

We introduce the maps a and q = (t, q) as in Lemma WA\ and we set E{s) := £(t(s), q{s)) 
for all s G [0,5"], so that E(t) = E{a(t)) for all t G [0,T]. Since t, g are Lipschitz 
continuous and £ is of class C\ the classical chain rule fl3.5p - fl3.7l) yields 

E'{s) > dtE,{t{s),q{s))?{s) - \dg8.\ {t{s),q{s))) \q'\{s) for £-a.a. s G (0,5). (4.18) 

On the other hand, since £■ is a Lipschitz map and since ^ = /Xq, the general chain rule 
of |AmD90j yields 

77^° = (^'oa)/i^° > (9i£(t,g(t))?oa- |9,£| (t,g(t)) 1^1 °^)/i™. (4.19) 

Taking into account (14. 9p . we conclude (notice that ii^' o a is well defined /iq°-a.e., since, 
for every Lebesgue negligible set A^ C Cq = cr(Cq), (14. 8 p yields ^'^{a~^{N)) = 0). ■ 

Now we are able to define the notion of BV solution. The formulation is more complicated 
than the one defining parametrized metric solutions, but it nicely reflects the different 
flow regimes of rate-independent flow, and the jumps. A shorter but much more implicit 
formulation will be given in Remark 14.41 



Definition 4.3 (BV solution) Let (Q,rf, £) satisfy (KQ, flQ]) . flglE]) . A function 



q G BV([0, T]; Q) is called a BV solution of the rate-independent system (Q, d, £), if the 
following four conditions hold: 



£(ti, g(ti)) - £(to, g(to))- / ' dtEit, q{t)) dt 

J to 

<-So(g,[to,ti])/orO<to<ti<T; 



(4.20a) 



|a,£|(t,g(t))<l /ortG[0,T]\J,; (4.20b) 

\d,E\ {t, q{t)) > 1 forte sppt(;u,); (4.20c) 

for t G Jg there exist y* G A{q(t^), q{t^)) and 6^ G [0, 1] such that 

(a) y\e') = qit), 

{(3) \d,8,\{t,y\e)) >1 for all ee [0,1], 

(7) £(t,g(t+))-£(t,g(r)) = -/,^|9,£|(t,y*(0)) \y'mde. 



(4.20d) 



Again we point out that, due to the chain rule inequality (I4.16p . relation (I4.20ap holds 
as an equality, which is the energy balance. Using this energy identity on the intervals 
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[t—h, t] and [t, t+h] and letting h\0 leads to the first two of the following jump relations, 
which will be used later (recall the definition fl4.1ip of 6*1): 

£(t,g(t-))-£(t,g(t)) = So{t;qit-),q{t)) = Si{t;q{t-),q{t)), 
£(t,g(t))-£(t,g(t+)) = So{t;q{t),q{t+)) = S,{t;q{t),q{t+)), (4.21) 

£(t, g(r ))-£(t, g(t+)) = So{t; q{t-), g(t+)) = 5i(t; g(r ), g(t+)), 

for each t G Jg. The third relation follows from (14.20d[) . By the definition of the slope 
distance 5*0, these jump relations already include the existence of a connecting gradient- 
flow curve y G A{q{t-) , q{t+)) , i.e. (a), {(3), and (7) of K2Mf follow. 

The above formulation of BV solutions looks quite lengthy compared to the more 
elegant forms of gradient-like flows, which can be characterized by one inequality, cf. e.g., 
(13. 8p or (I3.20cl) . However, this formulation reflects the mechanical interpretation of the 
three different flow types quite well, namely sticking, slipping and jumping. The following 
result presents a more compact form, which is however less tractable for further analysis. 



Proposition 4.4 In the setting of p.QP , fIS.Rp . fl3.Ep . let Si(-, [to^^i]) be the functional 
defined on BV([0, T], Q) via 

Ei(g, [to, ti]) := / ' max {|9,£| (t, q{t)) , l}/x^°(dt) + / ' (\d,8.\ (t, q{t)) - lY dt 
+ Si (to; g(to), q{tt)) + Siiti- q{t^), g(ti)) 

tGJqn(to,ii) 
Then, q G BV([0, T];Q) is a BV solution if and only if 

£(ti,g(ti))-£(to,g(to))- / ' 9t£(t,g(t))dt < -Si(g, [to,ti]) for < to < t^ <T . (4.22) 

Jto 



Proof: It is clear that under conditions fl4.20bp . (I4.20cp . and fl4.20dp a BV solution q 
satisfies 

So(g; [to, h]) = Si(g; [to, ti]) for < to < ti < T, (4.23) 

so that fH:^Uip yields (1123) • 

Conversely, if g G BV([0, T]; Q) satisfies (I4.22p . the chain rule (I4.16P and the inequality 
So(-; [to,h]) < Si(-; [to,ti]) yield (H:20al) and (14:23|) . Choosing e.g. to = 0,ti = T we get 

= j (max{|9,£| (t,g(t)) ,1} - \d,E,\ (t, g(t)) )^™(dt) + j {\d,E\ it,qit)) - l)^dt 
+ (S,{0; g(0), g(0+)) - 5o(0; g(0), g(0+))) + (s,{T; q{T-), q{T)) - So{T; q{T-), q{T)) 



+ 



Y, [Si{t; q{t-), g(t)) - 5o(t; g(t-), g(t))+5i(t; g(t), g(t+)) - 5o(t; g(t), g(t+))" 
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Since each addendum is nonnegative, we easily find ( 14.20b[) and f l4.20cp (recalling that 
|(9£| is continuous). Moreover, passing to the limit in (14.221) as to /" t, ti\ t, with t G Jg, 
we conclude 

£(t, g(t+)) - £(t, q{t-)) < -5i(t; g(r ), q{t)) - S,{t; q{t), g(t+)). 

Recalling (14.111) and the chain rule, for every t G J^ we find a curve y* satisfying condition 
gMM- ■ 

BV and parametrized metric solutions. We claim that the notion of BV solution 
is essentially the same as that of parametrized metric solution. Intuitively, (I4.20ap cor- 
responds to (13.231) . Further, ( ]4.20bl) and (I4.20cp are the analog in the BV setting of the 
first of (I3.24cp . which encompasses both sticking and rate- independent evolution (recall 
Remark 13. 5p . The jumping regime is accounted for by condition (]4.20dp : at jump times, 
the system switches to a viscous, rate-dependent behavior, following a path described by 
a generalized gradient flow, see (7) in (14.20dp . 

In order to formalize these considerations, we return to the trajectories in Q^. Indeed, 
we may associate with each BV solution gev a trajectory, by filling the jumps of the graph 
{ (t, QBvit)) I t G [0, T]) } with the curves y* G AC([0, 1], Q), for t G Jqbv Thus, we obtain 

T = { (t, gBv(i)) 1 1 G [0, T] } U |J { (t, y\9)) \9e[0, 1] }. 

iGJcjBV 

By construction, T is a connected curve that has exactly the length Var(g, [0,T]) + T 
if we use the extended metric dQr^{{to,qo),{ti,qi)) := |to— ti| + d{qo,qi) on Q^. Hence, 
there exists an absolutely continuous parametrization of T, and it can be shown that 
this parametrized curve is a parametrized metric solution. Indeed, in Example 17.41 we 
shall show, that to a given BV solution, there may correspond infinitely many distinct 
parametrized metric solutions. 

On the other hand, we can pass from parametrized metric solutions (t, g) defined in 
[so, si] to BV solutions by choosing 

a{t) G { s G [so, si] I t{s) =t} and defining q{t) := g((T(t)) . (4.24) 

Hence, Jg = { t G [0, T] | cr(t+) > cr(t^) }, and we see that q(t) is uniquely determined 
from g for t G [0, T] \ Jg. At the jump times t we can in fact choose any point q(t) = q{s) 
with s G [a{t-),a{t+)]. Note that 

y\9):=q{a{t-)+9[a{t+)-a{t-)]), 9e[0,l], t G J, (4.25) 

defines a connecting jump path as desired in fl4.20dp . We collect these remarks in the 
next proposition, whose proof easily follows from Lemma [4.11 (see also the Remark 13. 3p . 



Proposition 4.5 In the setting of I ^J^ , i KWtf . I KF^ . let ggy e BV([0,r]; Q) be a BV 



solution of the rate-independent system (Q, d, £) and let q = (t, q) be a map as in Lemma 



4-1 Then, setting 



z{s):={ '^'^ ^/^^C,,^ 

^ y\9) tfseJg, t{s)=t, s= {1- 9)a{t-) + 9a{t+) for 9 e [0,1], 
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the map {t,z) : [0,5*] ^ Qt is a parametrized metric solution of (Q,(i, £) according to 
Definition [XE 

Conversely, if {t,q) : [sq, Si] -^ Qt is a surjective parametrized metric solution (i.e. 
t(so) = and t{Si) = t), then any map q defined as in fl4.24p is a BV solution. 

The next result shows that BV solutions can be directly obtained as a vanishing 
viscosity limit, as in Theorem I3.8[ but now rescaling is not needed. The imposed a priori 
bound on the total variation for the viscosity solutions g^ can be easily obtained from the 
energy inequality (13.81) under general assumptions on £, see e.g. 



Corollary 4.6 (Vanishing viscosity limit (II)) Let{Q,d,E,) satisfy (|3.QD, (13. R|) . and 



(13. Ep . For every e > Q let q^ E AC([0,T]; Q) he a solution to (13.81) for ip = ip^- Assume 
that ge(0) ^ go ^s £ \ and Var(ge, [0,T]) < C for all e > with a constant C 
independent of e. Then, there exist a subsequence g^^. with e^ \ and a BV solution q 
for {Q,d, £) such that q£,^{t) — >• g(t) as k —>■ oo for all t G [0,T]. 

Proof: Let us consider the functions cr^ as in Lemma [4. 1[ By Helly's selection theorem 
we can find subsequences {qejk, {o'ek)k converging pointwise in [0,T]. Let us consider the 
corresponding parametrized metric solutions {te,q^e) introduced in Proposition 14. 5[ Since 
air is absolutely continuous with cr^ > 1, differentiating the identity cri;(t) = ^^(^^(t)) + 
Pe{(^e(t)) we obtain rrie := t'^ + |g^| = 1 a.e. in {0,3^) and S^ := (^^(T) > T. Since S^^, 
converges to 5* > T > 0, up to a further linear rescaling it is not restrictive to assume 
that Ssf. = S and me^. = Sef./S — > 1. 

Applying Theorem 13.81 we can find suitable subsequences (still labelled Ek) such that 
(tsk^'lek) -^ (^^) ^^ C°([0,S'];Qt)). Since ge(t) = g;(cre(t)) and t = fe((Je(t)),^e easily 
get ge^(t) -^ g(cr(t)) and t(cr(t)) = t, so that g is a BV solution induced by (t,q) as in 

mM- ■ 

We conclude the section with the Proof of Lemma I4.lt 
Part A) is immediate. Part B) is an obvious extension of [Fed69., 2.5.16], since each couple 
of points in Q can be connected by a geodesic. Notice that T{Vq{t)) = t and therefore 
f (p(cr(t))) = t if t G Cg. We thus get r o p = t in Cg. 

In order to prove C) it is not restrictive to assume A = [0,T] and ( G C°([0,T]). 
Then, (i2D follows from (|i3]) and |Fed69l 2.5.18(3)], observing that 

L{{se [0,S] \t{s) <t}) =Vq{t) foralltGCg, L{{r e [0,R] \?{r) <t}) =Vy{t). 

Let us now prove the first identity of (14. 9p : setting Jg := t^^(Jq) = [0, S] \ Cg, we observe 
that t (s) = for £-a.e. s G Jg. Since t is Lipschitz continuous and monotone, the change 
of variable formula and (14.70 yield, for every continuous function ( with compact support 
in(0,T), 

C(t)dt= / C{t{s))?is)ds= [ at{s))?is)ds= [ C(i)?oa(t)/i^°(dt). 

Jo JCa JCa 
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The second identity of (14.91) follows by a similar argument: 

C(t)^™(dt) = / at)f,,{dt) = [ C(r(r))dr = / C(r(p(s)))p'(s)ds 

JCq Jo 

The identity p' o a = |g'| o a follows from the property Vg(t) = Vg{a{t)) for all t G Cg, so 
that Vq{s) = p(s) for all s G C^. ■ 

5 Other solution concepts 

Here we discuss other notions of solutions for rate-independent systems (Q, d, £), namely 
energetic solutions, local and approximable solutions, and ^-minimal solutions. 

5.1 Energetic solutions 

The concept of energetic solutions provides the most general setting, in the sense that 
it does not even rely on a differentiability structure like the Finsler metric CR, but only 
uses the distance d. In such a framework it is even possible to consider quasi-metrics (i.e. 
unsymmetric and allowed to take the value oo), cf. |Mie05] . 

Definition 5.1 A mapping q : [0,T] ^ Q is called energetic solution for the rate- 
independent system {Q,d,8.) if for all t G [0,T] the global stability (S) and the energy 
balance (E) hold: 

(S) VgGQ: £(t,g(t))<£(t,g) + rf(g(t),g); 

(E) e(t, q{t)) + Var(g, [0, t]) = £(0, g(0)) + /J 9,£(s, q{s)) ds . 

We refer to |MiT99l IMTL02J for the origins of this theory and to |Mie05j for a survey. In 
analogy with (14.211) we have the jump relations 

£(t,g(r))-£(t,g(t)) = %(r),g(t)), 

£(t, g(t))-£(t, g(t+)) = d(g(t), g(t+)), (5.1) 

£(t,g(r))-£(t,g(t+)) = %(t-),g(t+)), 

for all t G Jg. Here they are easily obtained by considering the energy identity £(s, q{s)) + 
Var(g, [r, s\) = £(r, q{r)) + j\ dr^{T, g(r))dr, which follows immediately from (E), for the 
intervals [t—h,t], [t,t+h], and [t—h,t+h], respectively, and letting h\0. 

To compare energetic and BV solutions, we introduce the global slope S[£(^, ")] • ^ ^ 
[0, oo] via 

G\Pf+ \]( \ (£(t,g)-£(t,g))+ 11 ^ n 

g £(t, ■) (g) := sup^ — xr — for all gG-D. 

qH d{q, q) 
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Using this definition, the global stability (S) can obviously be rephrased as 9[8,(t, ■)]{qit)) < 
1. We also have 

|9,£|(t,g) < g[£(t, ■)](?) for all {t,q) G Qt • (5.2) 

Indeed, choosing a local coordinate system in Q, one can check (cf. [BCSOOt Ch.VI.2]) 
that 

ia,£|,,,). .up (M(^.u,„.,p(^Mz4M)):. ,5.3) 

^eT,Q\{0} ^i{Q,v) q^q d{q,q) 

whence (I5.2p . 

Remark 5.2 It is well known that (S) implies the lower energy estimate 
£(s,g(s))-£(r,g(r))- [ 9t£(t, g(t))dt > - Var(g, [r, s]) 



for < r < s < T, cf. |MTLn2[ Thm.2.5] and |Mien5l Prop. 5.7]. In the present setting 
this is in fact an easy consequence of the chain rule inequality (I4.16p and of the observation 
that g[£(t, ■)]{qit)) < C implies r(g, [r, s]) < C Var(g, [r, s]). 

Moreover, it is possible to derive a gradient-flow like inequality of the type given in 
(El]), dS^Qc]), or f02|) . For this, define the functional r,(-, [r, s]) on BV([0,T], Q) via 



r,(g,[r,s]):= / max {g[£(t, ■)](g(t)), lK°(dt) + / {9[Eit,-)Mt)) - if dt 

J r J r 

+ max{|£(r, g(r)) — £(r, g(r^))|, d{q{r), q{r^))} 
+ max{|£(s, g(s"))-£(s, g(s))|, rf(g(s"), g(s))} 

+ ^ [max{|£(t,g(r))-£(t,g(t))|,d(g(t-),g(t))} 

te(r,s)nJg 

+ max{|£(t, g(t))-£(t, g(t+))U(g(t), g(t+))}] . 

Then, g : [0, T] ^ Q is an energetic solution for (Q, d, £) if and only if S[£(0, •)](g(0)) < 1 
and 



^{s,q{s)) - 8.{r,q{r)) 



/ 9i£(t,g(t))dt< -r,(g, [r,s]) for 0<r < s<T. (5.4) 

J r 



The following result essentially states that every energetic solution g is a BV solution 
outside its jump set. Moreover, if the jump relations (14.211) and (15.11) are both satisfied, 
then an energetic solution is also a BV solution. Conversely, if a BV solution additionally 
satisfies S[£(t, ■)]{q{t)) < 1, then it is an energetic solution as well. 

Proposition 5.3 (Comparison between energetic and BV solutions) Assume that 
(Q,c/,£) satisfies I ^J^ , (lOD . and (ISlEl) . 

(A) If q & BV([0, T]; Q) is an energetic solution, then q is also a BV solution if and only 
if ( I4.20d|) holds additionally. 

(B) //gGBV([0,T];Q) is a BV solution with 9[Z{t,-)]{q{t)) < max{l, |9g£| (t, g(t))} for 
all t G [0, T] and S[£(0, ■)]{q{0)) < 1, then q is also an energetic solution. 
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Proof: To prove (A), we first note that the necessity of ( I4.20d[) is obvious at it is part of 
the definition of BV solutions. To establish the sufficiency in (A), we observe that (]4.20dp 
yields (14.211) . so that the dissipation term Si of Proposition 14.41 satisfies 

Si(g, [to,h]) < r,(g, [to,ti]) for all < to < ti < T. 

Hence, (14.221) follows from (15.41) . Thus, statement (A) is established. 

The necessity of the additional condition on S[£(t, ■)]{<l{t)) is obvious, since energetic 
solutions have to satisfy the stronger stability condition (S). To show the sufficiency we 
observe that the additional condition yields 

r.(g, [to,h]) < Si(g, [to,ti]) for all 0<to<h<T, 

so that ([EID follows from IK2T\ . ■ 

Remark 5.4 The additional condition in Proposition 15.3( B) is implied by the general 
condition 

g[£(t, •)](?) = \d,e\ (t,g) for all (t,g) G Qt ■ (5.5) 

If this condition holds, then the notions of energetic solutions and BV solutions coincide 
under the additional assumption that the initial state go is stable, i.e. S[£(0, ■)]{qo) ^ 1- 
One condition guaranteeing (15.51) is a metric version of convexity for £(t, •), see |AGS05i 
Def. 2.4.3]. Here, we say that 5" : Q ^ M U {oo} is convex on (Q, d), if 

Vgo,gieQ, 9 e [0,1] ^qe e Q : d{qo,qe) = Od{qo,qi), d{q0,qi) = {l-e)d{qo,qi), 

3'{qe)<{l-e)^{qo) + e3^{q,). 

To establish (15. 5p for 3^ note that for each g and e > we have S[?'](g) > 2^~) ~g 
for some q. Moreover, for each n E N, choosing 9 = 1/n we find g„ with d{q, qn) = -d{q, q) 
and d{qn,q) = '^^d{q,q). Applying (15.61) we obtain 

\Oq-^\{Q) > hmsup — —z — - — > hmsup " 

n^oD d[q, qn) n^oo — d[q, q) 

J{q)-J{q) 



d{q,q) 



>S[3-](g)-e. 



In this way, part (B) of Proposition 15 . 31 is a generalization to the metric setting of Theorem 
3.5 in |MiT04j . which states that for a Banach space Q, a convex energy functional £, 
and a translation invariant metric d the sub differential formulation and the energetic 
formulation for (Q, d, £) are equivalent. 

5.2 Local and approximable solutions 

As we have already mentioned in the introduction, energetic solutions have the disadvan- 
tage that the stability condition (S) is global, so that solutions tend to jump earlier than 
expected, see Example 17.11 To avoid these early jumps, the vanishing viscosity method 
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was employed in |EfM06l IDD*07t IToZOGt IKMZ07J . When avoiding parametrization and 
studying the hmits of the viscous solutions q^ : [0, T] ^ Q directly, one obtains an energy 
inequality and a local stability condition. Hence, we next introduce the notions of local 
solution and of approximable solution, generalizing the definitions given in |ToZ06j to the 
metric setting. 

Definition 5.5 A mapping q : [0,T] ^ Q is called local solution, i/ (al) and (a2) hold: 

(al) \dge\ (t,g(t)) < 1 for a.a. t e [0,T]; 
(a2) for all r,s & [0, T] with r < s we have 

£(s, g(s)) + Var(g, [r, s]) < £(r, g(r)) + // 9t£(r, g(r)) dr. 

We will see in the examples of Section [7| that the notion of local solution is very general. 
Using (15.21) . it is clear that all energetic solutions are local solutions. Similarly, all BV 
solutions are local solutions. To see this, we use fl4.20cp and (]4.20dl) to obtain (al), since 
there are at most a countable number of jump points, and we conclude So(g, [r,s]) > 
Var(g, [r, s]) for < r < s < T, which gives (a2). 

On the other hand, note that, unlike the case of energetic solutions, the combination of 
the local stability condition with the energy inequality does not provide full information 
on the solution q. In particular, the behavior of the solution at jumps is poorly described 
by relations (al) and (a2). This also highlights the role of the term So(g, ■), here missing, 
in the energy identity for BV solutions. As a consequence there are many more local 
solutions, see also Example 1 7. 1[ 

Indeed, the vanishing viscosity method turns out to provide a selection criterion for 
local solutions. Among local solutions, we thus distinguish the following ones: 

Definition 5.6 A mapping q : [0,T] — > Q is called approximable solution, if there exists 
a sequence {ek)k&n with Sk \ and solutions g^^. G AC([0,T], Q) of (13.81) with ip = ipsk 
such that for all t G [0, T] we have qe^if) — > q{t). 



It follows from Corollary 14.61 that, under the assumptions ( 3.Q ), fl3.R|l . and fl3.Ep . any 
approximable solution is a BV solution as well. 

The notion of approximable solutions suffers from two drawbacks. First of all, there is 
no direct characterization of the limits in terms of a sub differential inclusion or variational 
inequahty, unlike for parametrized/BV solutions, recall Proposition 13. 7[ Secondly, since 
the solution set is defined through a limit procedure, it is not upper semicontinuous with 
respect to small perturbations, as shown in Example 17.31 This is in contrast with the 
stability properties of the set of parametrized/BV solutions, see Remark I3.1UI 

5.3 Visintin's ^-minimal solutions 

In |Vis01t IVisOG] a new minimality principle was introduced. Here, we present the adap- 
tation to rate-independent evolutions proposed in |Mie05t Sect. 5. 4], in the current metric 
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setting. Again, we use parametrized curves, as it is essential to have continuous paths. 
For simphcity, we restrict to arclength parametrization, i.e., 

t'{s) + \q'\{s) = 1 for a. a. s G [so,Si]. (5.7) 

In the framework of p.QP , (13. Rp . (13. Ep . for a given initial pair (to,qo) we introduce the 
space of arclength-parametrized paths (on some interval [0, S]) starting in (to, Qo) via 

■^^(to, go) := { {t, q) e C°([0, S], Qt) I t(0) = to, t nondecreasing, 

g(0) = go, t{s) + Var(g, [0, s]) = s for all s G [0, S] }. 

On this set we define the function $ : Asito, go) -^ L°°([0, S]) via 

$[t,g](s) = e(t(s),g(s)) + Var(g,[0,s])- / 9i£(t(r), g(r))t'(r)dr . 

Jo 

Between paths in ^^^(to, go) we introduce an order relation ^ as follows. For (t, g), (r,p) G 
-^5(^0, go), define the "arclength of equality" via 

S[(t,g),(r,p)] = inf{s G [0,5] | (t(.),g(s)) ^ (r(.),p(.)) }. 

Then, the order relation is given by 

.^ .^. . ^^ f Vs>S[(t,g),(r,p)] 3 s, G(S[(t,g),(r,p)],.): 

1 $[(t,g)](s.)<$[(r,p)](s.). 

Definition 5.7 74n arclength-parametrized function {t,q) : [so, Si] ^ Qy is called a $- 
minimal solution /or (Q, (i, £), if for all {t,p) & As^-SQ{t{sQ),q{so)) we have (t,q) :< {t,p). 

Like for energetic solutions, this solution notion appears particularly suitable to handle 
nonsmooth energy functionals, since no derivatives/slopes of £ with respect to the variable 
g occur in the definition of the functional $. 

We now show that, in a smooth setting, ^-minimal solutions are parametrized metric 
solutions. Using suitably chosen test functions, it can be shown that a necessary condition 
for ^-minimality is the local condition 

— $[(t, g)](s) < K(t(s), q{s)) for a.a. s G (sq, Si), 

where X : Qy — > K is defined via 

y{{t, q) = liminf (- inf{ £(t, g)-£(t, q)+d{q, q) \ d{q, q)<e}]. 

A simple calculation gives 

m,q)^i ' for|.,£|(t,g)<l, 

^ l-|9,£|(t,g) for|5,£|(t,g)>l 
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Since t' + |g'| = 1 a.e. and 

-^$[(t, q)]is) = ^ms), qis)) + \q\is) - 9,£(t(s), g(s))t'(s) for a.a. s G (so, Si), 
els as 

we conclude that all <l>-niinimal solutions satisfy a.e. in (sq, si) 

^£(t(s),g(s))-a,£(t(s),g(s))t'(s)<-M(t'(s),|g'|(s),|9,£|(t(s),g(s))), 

together with the constraint t'{s) + \q'\{s) = 1, where M{a, i^, = ^ + (^ ~ l)"*"- We have 
thus proved that any ^-minimal solution (t, q) on [sq, Si] is a parametrized metric solution, 
and hence a BV solution (up to a parametrization). The opposite is in general not true, 
see Example 17.21 Further, Example 17.31 shows that the set of ^-minimal solutions is not 
stable with respect to perturbations. 

6 Outlook to the analysis in metric spaces 

In |MRS08j we shall analyze rate-independent evolutions in 

a complete metric space (X, d) . (6.1) 

In fact, using the results in |RMS08] we shall be able to handle the case in which d is a 
quasi- distance on X, i.e. possibly unsymmetric and possibly taking the value cx). 

In this framework, the metric velocity (11.31) of a curve q G AC([0,T]; X) is defined by 

\q'\{t) ■.= \im -d{q{t),q{t+h)) = lira -d{q{t-h),q{t)) for a.a. tG(0,T). (6.2) 

h\0 h h\0 h 



In the Finsler setting (|3.QD-(I3.RP of Section [3TT| one indeed has \q'\{t) = 3^i(g(t), g'(t)) 



for a.a. t G (0,T) (see |BCS00l Chap. VI. 2]). Further, given a functional 

£ : [0, T] X X ^ M U {oo}, with domain dom(£) = [0, T] x D (6.3) 

for some D G X, the local slope is defined via 

|9,£|(t,g)=limsup-i-(£(t,g)-£(t,g))+, (6.4) 

q-^q d{q, q) 

which in the Finsler setting X = Q coincides with (13.41) . With these tools, (13.81) is the 
purely metric formulation of doubly nonlinear equations of the type (13. 9p . provided the 
local slope fulfils the following chain rule inequality. 

Definition 6.1 We say that the triple (X, (i, £) satisfies the chain rule inequality if for 
every absolutely continuous curve (t,q) : [sq, Si] — > [0,T] x X such that t' > a.e. in 
(so, Si), g(s) G D for every s G [sq, Si], and 

' |9,£| (t(s),g(s)) |g'|(s) + |9i£(t(s), g(s))| t'(s)ds < oo , 

'so 

the map s ^—>- £(t(s),g(s)) is absolutely continuous on [sq, si] and satisfies 
d^ 



£(t(s),g(s))-9i£(t(s),g(s))t'(s)>-|9,£|(t(s),g(s))|g'|(s) a.e. zn(so,si). (6.5) 
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Unlike in Section [3?ll within this abstract setting the chain rule inequality is no longer 
granted, but has to be imposed instead. We refer to |RMS08j for a discussion on some 
sufficient conditions for (16.51) to hold. 

If the above chain rule holds, the (parametrized) metric formulation (I3.12p is again the 
starting point for the vanishing viscosity analysis, which was developed in Section 13.21 in a 
Finsler setting for smooth £. In this general framework, £ has to satisfy some coercivity 
and (lower semi-) continuity properties: 

the functionals £(t, ■) are lower semicontinuous and 

uniformly bounded from below with Kn := inf £(t, o) > — oo : 

Vt G [0,T] : £(t, ■) has compact sublevels in X; 

3Ki>0VgGZ^: £(-, g) G C^([0,T]) and (6-6) 

|ai£(t,g)| <iri(£(t,g)+l) for alH G [0,T]; 

V((t„,g„))„,6N C [0,r] X X with (t„,g„) -^ {t,q) : 

dt8,(t,q) = lim dt8.{tn,q„) and \dg8.\ (t,q) < liminf |9g£| (t„,g„) . 

In |MRS08j . under assumptions (16.11) and (I6.3l) - (l6.6p on £, we shall perform the vanishing 
viscosity analysis of Theorem 13. 8^ leading to the notion of parametrized metric solution of 
the rate-independent system (X, d, £). In this general setting, it is obviously still possible 
to consider the notion of BV solution, and our remarks on the comparison between BV 
(parametrized) and local/approximable/$-minimal solutions carry over. 

Indeed, in |MRS08j we shall discuss BV solutions with more detail, in particular prov- 
ing existence through approximation by time discretization and solution of incremental 
(local) minimization problems. 

7 Examples 

Many of the differences between the various solution concepts discussed above manifest 
themselves already in the case in which the state space is the real line. Hence, we discuss 
the very simple model with 

Q = M, d{qo,qi) = \qo-qi\, e{t,q) = U{q) - i{t)q , (7.1) 

where the function i will be specified in the different examples. The potential U is the 
nonconvex function given via 



U{q) = { A-W for |g| < 2, (7.2) 



\{q+^? 


for g < -2 


^-w 


for q\ < 2, 


W-^f 


for g > 2. 



As initial datum we shall take 



go = -5. (7.3) 
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Example 7.1 We let i{t) = t for all t > 0. We claim that the approximable, the $- 
minimal, the parametrized, and the BV solutions on [0, cx)) are essentially unique and 
coincide. However, the unique energetic solution is different. Moreover, we show that 
there is an uncountable family of different local solutions. With direct calculations, one 
sees that the energetic solution takes the form 

q{t) = t-5 for t e [0, 1) and q{t) = t+3 for t > 1. 

Choose any t^, G [1,3] and any q^ G [3+t*, 3+t*+min{2,4A/t^, — 1}] . Then, 

{t-5 fort G [0,t), 
q^ fort G (t,g*-3], 
t+3 for t > g*— 3, 

is a local solution. Note that the starting point of the jump at q{t^—) = t*— 5 can be 
chosen in a full interval. Moreover, for a fixed t^: > 1 we still have the possibility to choose 
the ending point q^, = q{t^:+) of the jump in a full interval. 

All the other solution types essentially lead (up to definition in one point) to the same 
solution. Without time parametrization it reads 

t-5 forte [0,3), 
lit) = { q^ fort = 3, 
t+3 for t > 3, 

where g* G [—2, 6] is arbitrary. The associated arclength-parametrized solution takes the 
form 

( (f,f-5) for sG [0,6], 

{t{s),q{s)) = l (3,s-8) for sG [6, 14], 

[ (i-4,f-l) for. > 14. 

Example 7.2 In this example we show that, in general, approximable solutions and $- 
minimal solutions are different. In particular, recalling the discussions in Sections \5.2\ and 
\5.3\ this shows that the set of BV solutions (or parametrized metric solutions) is strictly 
bigger then any of the other solution sets. 

In the setting of fl7.ip - fl7.3l) . we now choose the function i{t) := min{t, 6— t} for all 
t > 0, i.e., the loading reduces exactly when the solution reaches the jump point. It is 
easy to see that there are two different BV solutions: qi, which jumps at t = 3, and q2, 
which does not jump. We have 



Qiit) 



t-5 for te [0,3), 

6 for te (3,5], 

11-t for te [5,9), 

3-t for t > 9; 



t-5 forte [0,3], 
q2{t) = { -2 forte [3,5], 
3-t for t > 5. 
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For e > the viscous solution (f of the differential inclusion 

G Sign(g) +eq + U'{q) - £(t), g(0) = -5, 
is unique and can be calculated by matching solutions of linear ODEs. We find 

{t-5+e(e-*/=-l) for te [0,3], 
ql fortG[3,t^], 

3-t+e(e-(*~**)/^-l) for t>tl, 

where ql = ^^(3— ) < —2 and tl = 3~ q^ > 5. Thus, we have q'^it) — > q2{t) for every t > 
as e I 0, and q2 turns out to be approximable, whereas gi is not. As a general principle, 
one may conjecture that viscosity slows down solutions, and thus approximable solutions 
tend to avoid jumps if there is a choice. 

For ^-minimal solutions this seems to be opposite. We claim that q\ is (up to a 
reparametrization) ^-minimal but q2 is not. For this, we use the arclength parametriza- 
tions 

r (f,f-5) for sG [0,6], ( (f , f-5) for sG [0,6], 

(ti,qi){s) = I (3,s-8) for s G [6,14], (t2,g2)(s) = < (s-3, -2) for s G [6,8], 

[ (s-11,6) for sG [14,16]; [ (|+1,2-|) for s > 8. 

The functionals ^j{s) = ^[(tj,qj)]{s) for all s >0, j = 1,2, can be calculated explicitly: 
indeed, one checks that 

ip,{s) = l } for sG [0,6], <^2(s) = -fors>0, (7.4) 

^^' \\-\{s-Qf for sG [6, 10], ^'^^2 - ^ ^ > 

which clearly shows that (^2, qi) Is not ^-minimal for s G [0, 7]. 

To prove ^-minimality of {ti,qi) we point out that chain rule inequality ( \3.7\i gives 

HT{s),pis)) > ^ + Var(p, [0, s]) - ^ |9/| (r(a),p(a)) \p'\ia) da (7.5) 

for all s G [0, T] and all {t,p) G At{0, go)- Equality holds in (17. 5p if and only if {t,p) is a 
parametrized metric solution (t,q) on [0,T]. In that case, in view of (I3.24al) one further 
has, for all s E [0, T], 

Ht{s),q{s)) = ^ + Var(g, [0,.]) - j^ M{\?\{a), |g'|(a), \d,8,\ {t{a),q{a)) WHa) da. 

Therefore, in order to check that {ti,qi) is ^-minimal, it is sufficient to prove that 
(ti,Qi) ^ (tjO) ^or all parametrized metric solutions {t,q), and this, for all the arclength 
parametrizations {t,q) corresponding to the (not jumping) BV solution q2, follows from 
the previous discussion on (t2,q2). Now, the above energy balance states a general fact 
about parametrized metric solutions: $(t,g) is constant as long as no jumps occur, i.e. 
\dq8.\ (t,q) < 1 holds. If jumps with \dgE,\ (t,q) > 1 occur, then $ will strictly decrease. 
Thus, if there is a choice between one solution with a fast jump and another without 
jumps, then the solution without jumps cannot be ^-minimal. 
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Example 7.3 Here, we study the parameter dependence of solutions under the loading 

£s{t) = mm{t, Q+26-t} for t > 0, 

where 6 is a small parameter. In the case 6 = we have two BV solutions qi and q2 (or 
similarly parametrized metric solutions), as was discussed in Example \7.2[ For — 1 < 5 < 
there is only one solution, namely 

t-5 forte [0,3+6], 
q'it) = <j 5-2 forte [3+6,5+5], 
3+26-t for t > 5+6. 

The corresponding parametrized solution is the unique ^-minimal solution. Now, for 
6 y we find q^{t) — > q2{t) for every t > 0. Hence, the set of ^-minimal solutions is 
not closed (or ''not stable" or "not upper semicontinuous" ) under pointwise convergence. 
Similarly, we may consider 6 > to obtain a unique BV solution q^ that jumps at 
time t = 3 before the unloading starts at t = 3+6 > 3. Clearly, these solutions are 
approximable and converge pointwise to qi, which is not approximable. Thus, the set of 
approximable solutions is not upper semicontinuous. 

Example 7.4 We provide an example where one BV solution corresponds to many dif- 
ferent parametrized metric solutions. The BV solution has exactly one jump, and there 
are infinitely many distinct connecting orbits y in (iv) of Definition \4.3[ giving rise to 
infinitely many distinct parametrized metric solutions. We consider 

Q = R^ and d{q,q) = -{[qi -qi[ + [q2 -q2[) ■ 
With q = {qi,q2) e Q = M^ the potential takes the form 

where U is defined in ( Fr2l) and W : R ^ [0,oo) by W{p) = for [p[ < 1 and W{p) = 
(IpI — 1)^ else. Starting from g(0) = (—5, —5), we have q(t) = {q(t), q{t)), q being the BV 
solution of Example \7.1[ Hence, the (unique) jump occurs at t = 3, starting in (—2, —2) 
and ending in (6, 6). However, the set of connecting paths y is infinite. Indeed, for every 
connecting path there holds for a. a. s e {sq, Si) 



jj, I yi[s) + y2{s) 



ll/'l(^) = ^(ll/'i(^)l + 11/2(^)1), m{t,-)[{q{s)) = \ 

if |yi(s)— 1/2 (s) I < 1. Now, for a given curve 'j : [0, 1] -^ [0, oo) let us sety^ := (g— 7, q+'j)- 
Indeed, [y'^[ = l/2{[q' — ■y'[ + [q' + 7'|) = [q'[ whenever |7'| < [q'[. Therefore, 

|9,£(t, ■)[iy,is))[y'^[is) ds = f [d,E{t, ■)[{y{s))[y'[{s)ds 

for all curves 7 with 7(0) = 7(1) = and |7'|(s) < Ig'l(s) for a.a. s e (sq, Si), and for such 
•y's y^ is an optimal connecting curve. 
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